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We show that a set of N identical qubits coupled to a single cavity resonator can be decomposed 
into a set of independent subsystems, which can be thought of as higher pseudospin generalisations 
of the Jaynes-Cummings model. We derive and analyse the solutions to the equations of motion 
and demonstrate unusual beating behaviour resulting from a new form of -yn-type non-linearity 
appearing within a pseudospin ladder. Furthermore, we propose a relative phase shift transformation 
which allows one to switch/initiate the multi-qubit state in a desired pseudospin configuration, and 
show how such transformations can be used to undo spontaneous single qubit decay. We discuss the 
conditions which justify the validity of the rotating wave approximation in an iV-qubit system. 
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Superposition and entanglement are among the most 
striking characteristics that distinguish quantum from 
classical physics and give way for quantum information 
processing (QIP) [1, 2]. The building block behind the 
QIP technologies is a two level system (a qubit), and 
quantum algorithms rely on coherent manipulation of a 
single qubit and on controllable coupling between qubits. 
Recent years brought a tremendous progress in manipu- 
lation and coherence control of multiple qubits [3,4]. The 
control of the individual qubit as well as its read-out is 
often done by means of coupling it to a resonator (an op- 
tical [5] or a microwave [6, 7] cavity), in the framework of 
a cavity QED or a circuit QED architecture. In order to 
use the benefits of the speed-up offered by quantum pro- 
tocols, one needs to be able to involve a larger number of 
qubits and place them in a quantum register, to be able 
to control the coherence as well as the evolution of this 
multipartite system. Very recently, multiple-qubit regis- 
ters with three qubits coupled to a single resonator have 
been fabricated and coherent single-excitation exchange 
between them was demonstrated [8]. 

Interaction of a single qubit with a resonator is de- 
scribed by means of the Jaynes-Cummings model [9]. 
The behavior of this model is well understood. It pre- 
dicts such phenomena as the Rabi splitting and Rabi os- 
cillations facilitated by a change in the resonators popu- 
lation, quantum beats, and electromagnctically induced 
transparency [10]. 

The multipartite analog of the Jaynes-Cummings 
model, essential for the description of several qubits 
coupled to the same resonator, is known as the Tavis- 
Cummings model (TCM) [11] and is considerably less 
understood. The eigenvalues of the Tavis-Cummings 
Hamiltonian have been found [11], and the full dynamics 
of the two qubit system has been studied [12]. The in- 
homogeneous coupling has been addressed [13]. The un- 
derstanding of the model beyond two qubits is very lim- 
ited [14, 15], whereas the multitude of excitation schemes 



available in the model suggests rich and interesting dy- 
namics. 

In this Letter, we study TCM for N identical qubits 
resonantly coupled to a single resonator (bosonic mode). 
This system is symmetric under any qubit permutations, 
giving rise to a preferred state description and a system 
decomposition in terms of higher pseudospin composite 
particles (highly entangled multi-qubit states) . This per- 
ception of the problem allows to isolate the dynamics into 
a finite number of clusters of states and greatly simplifies 
calculations, reducing the number of dimensions from 2 N 
down to N + 1. We show that every unique pseudospin 
particle (a multiplet) is characterized by its own Rabi 
splitting and that single-qubit Rabi oscillations are re- 
placed by beating-like oscillation patterns. Only at the 
large number of excess resonator excitations these lead to 
the equidistant splittings observed in NMR experiments 
on spin multiplcts. Wc demonstrate furthermore that 
with a single multiple-qubit register one has access to the 
entire range of individual multiplets by means of condi- 
tional phase transformations, and show how this setup 
can be used to detect and fix single qubit decays (spin 
flips). 

The system. The TCM Hamiltonian reads 



H — H + Hi , 



N 
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Here, the operators a' , a refer to the resonator, and di 
describe the qubits. Assuming resonant u> = fijVi and 
equal coupling, gi = g, wc define the interaction picture 
Hamiltonian 



V = e iHot H ie - iHot = g (j+a + J_a f ) 
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The propagator U (t,f), describing the dynamics of 
the system under the assumption it is closed (very stable 
qubit and a non-dissipative cavity), is obtained by mul- 
tiple integration of the Schrodinger equation id t \tp (t)) = 
V\ijj(t)). For resonant coupling, this yields in the inter- 



very important. Here 



action picture U (t, 0) = e 
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evaluation of this sum is difficult, as for N qubits the 
number of dimensions to consider grows exponentially 
n<iim = 2 W . It would also obscure interesting physics of 
higher pseudospin particles. 

The pseudospin is formed in the following process. 
A set of N identical qubits has an iV-fold permutation 
symmetry: any pair-wise exchange of two qubits does 
not affect the system dynamics. The optimal (sym- 
metry preserving) basis can be determined from the 
irreducible representations (irreps) of the permutation 
group. The simplest example is that of two qubits 
(N = 2) spanned by computational basis |00), |01), 1 10) 
and 1 11). These can be combined into a spin singlet 
\j = 0,m = 0) = (|01) - 1 10)) /y/2 (antisymmetric) and 
a spin triplet \j = l,m=l) = |11), \j = 1, m = 0) = 
(|01) + 1 10)) /V2, \j = 1, m = -1) = 1 00) (symmetric un- 
der the exchange of the qubits), where we use the nota- 
tion adapted from the total angular momentum j and 
the magnetic m quantum numbers. For N > 2 one uses 
the Young tableaux to construct these states, see the 
supplementary material. In this formalism, the states 
\j = 0, m = 0), and (\j = l,m = or ± 1)) come from 
the ( | 1 1 2 1 ) diagram. 

Three qubits. The possible diagrams arc |i|2|3| , [Jp-l 
and corresponding to a quadruplet and two dou- 

blet states, reducing a single problem from 2 3 dimensions 
down to three problems in 4, 2, and 2 dimensions. Note 
that the states obtained from these irreps are not mutu- 
ally orthogonal, and need further orthogonalisation. 

In the supplementary material we present an algorithm 
for calculating the propagator for arbitrary initial condi- 
tions and initial qubits-resonator states correlations. Us- 
ing this protocol we calculate the three-qubit propagator 
in the basis spanned by a quadruplet \j = |) and two 
doublets \j = i). The solution to the doublet problem 
is a subject of standard quantum optics treatments [10], 
and the propagator for the j ' — | subspace in the basis 
\j,m) = {||, |), ||, |), ||, — |), | j,— §)} reads 



/ F ,i(n) F ltl (n)a F 2A {n)a 2 F 34 (n)a 3 \ 

a)Fi A (n) F 0:2 (n) F li2 (n) a F 2a (n) a 2 

{a)) 2 F 2 . l {n) tfF li2 (h) F 0|3 (n) Fi |3 (n) a 

V (at) 3 F 3 ,!(n) (at) 2 f 2 , 2 (n) &F lt3 (n) F 0A (n) J 



where the position (before or after) of the creation and 
annihilation operators relative to the functions F^j is 



2n+ 7 
A 



+ ,2 



cos(w±,2£) , 



and the rest of these functions as well as the frequencies 
w± a and the coefficients A± ; 2 are given by Eq. (7) in 
the supplementary material. The populations of different 
states are shown in Fig. 1. 

General case. For ./V qubits one can always block di- 
agonalisc the Hamiltonian, independent of the number of 
bosonic excitations or the coupling strength. As a result, 
we decompose this system with 2 N available states into 
a set of n C^n_ j , where [^J denotes the floor-function, 
i.e. one rounding down to the nearest integer. We label 
them with j being a (odd half-) consecutive integer val- 
ued for an (odd) even number of qubits involved, such 
that y > j > (i). Moreover, some j labelled multi- 
plets will be abundant more than others, with the their 
abundance Abjvj given by 
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The first set of block diagonalising states for two, three 
and four qubits, are listed in the supplementary material. 

Using this form we can predict how many blocks of 
each type will be present and calculate the within-the- 
block transition frequencies. The size of the largest 
(N + 1) x (N + 1) block grows linearly, as compared with 
the exponential growth in the computational basis, as the 
number of qubits increases. 

The 7V-qubit TC Hamiltonian decomposes into a num- 
ber of high pseudospin Jayncs-Cummings-likc terms in 
the j representation 

Mi} 

Hj c = SJj^ + g^^+a + Jj-tf) , 
j 

Jj,+ = X! v j>m \m+l){m\, J 7 -_ = Jj j+ , 

m=-j 

Jj,z = ^[Jj,+,Jj,-], Vj,m = VU~ m ) (1+j+m), 

where 5 is the qubit detuning from the resonator and 
<;, < q < N, is even (half odd) integer valued for an 

even (odd) number of qubits, e.g. Hj c = ( H^ c 
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For simplicity, we choose resonant coupling, 8 = 0. 
This reduces the problem of an arbitrary number of iden- 
tical qubits coupled to a single resonator to the calcula- 
tion of the propagator for any j-reprcsentation. Since 
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FIG. 1: (Color online). Population of the m = — 1,4) (dashed orange) 



■5,3) (dot-dashed red), J m = |, 2) (dotted 



blue), and \m = §, l) (continuous purple) states as a function of time (in units inverse coupling strength), given that the state 



is initiated in the m 



\, n = 4) state. 



j'-valuc H Degenerate States \j, n) 
1 |i, n >,|-a,n + l> 

1 |l,n), |0,n + 1), |-l,n + 2) 
|,n),|i,n + l), 
-i,n + 2),|-|,n + 3) 

2,n), |l,n + 1), |0,n + 2), 
-l,n + 3>,|-2,n + 4> 
f,n),||,n + l),|i,n + 2), 
_|,„ + 3>,|-|,n + 4),|-|, 



5) 



Spl itting [3] 



0,±2,/n- 



Large n Splitting [<?] 
sa ±v^ 
» 0, ±2vn 



±,/5(n±2)±4,/(n + 2) 



2 , 9 
4 



0,±yi0 (n+ f) ±6y^+ I) 2 + 2 ^ 0, ±2^, ±4^ 



n = Splitting [3] 
±1 

0, ±2.449 
±1.207, ±4.306 

0, ±2.787, ±6.499 



No analytical solution sa ±v / ™> ±3-v/", ±5^ ±1.364, ±4.744, ±8.979 

TABLE I: Summary of Rabi splittings. 



the Hamiltonian is time independent, the propagator is 
a (2j + 1) x (2j + 1) matrix taking the simple form of 

Uj (t) = exp (itHj C ^J . This can be straightforwardly cal- 
culated by diagonalisation, exponentiation and inverting 
the diagonalisation. The eigenvalues of Hj are the new 
Rabi field splittings, however they no longer correspond 
to the transition frequencies of the populations between 
the \j, m) and \j, m!) states, as we saw for propagators 
for j = i or j = 1 versus j = | pseudospin. 

We write the state of the j th pseu- 
dospin interacting with the cavity as \ipj) = 

Min(2j,n) 

2\2 Xi\ji m = — j ^ + i) \ n ~ i) y where the mini- 

i=0 

mum is taken as an upper limit of the sum since for a 
state ladder composed of 2j + 1 states, there might be 
fewer than 2j photons to climb the ladder. Thus, the di- 
mensionality of the state is restricted by the total initial 
number of excitations, i.e. if one initialises the system 
in the state \j = 4, m = — 3)|2), then the only reach- 
able states are \j = 4, m = — 4)|3), \j = 4, m = — 3)|2), 
\j = 4,m = — 2) 1 1) , and |j = 4, m = — 1) |0) . This means 
that both the pseudospin and the number of available 
excitations quantitatively and qualitatively affect tran- 
sitions frequencies between the available states. A list of 
eigenvalues of the hybridised states are given in Table 
I for 2j < n. In Fig. 2 we show for a nonet (j = 4 



state) how different splittings emerge as the number 
of available excitations grows. For 2j <C n {excess 
excitations, i.e. what is left upon climbing from the base 
to the top of the ladder) , all of the spacings of different 
pseudospins asymptotically converge to ±kgy/n, where 
k denotes an odd (even) integer, for an odd (even)-half 
integer valued j. In this semi-classical (large n) regime 
the shifted frequencies form a overtones-like spectrum, 
however in the regime of low n, they acquire \/n-like 
non-linearities which differ qualitatively for different 
pseudospin multiplets and also for incompletely used 
multiplcts, i.e. for 2j > n. These properties are pictured 
in Fig. 6 in the supplementary material. 

Switching pseudospins. The permutation symmetry 
breaks the complete set of states into submanifolds of dif- 
ferent pseudospins, which cannot be switched by means 
of the evolution of the system alone. There exists how- 
ever a multiply-conditioned phase-shift transformation 
which allows one to manually change between the sub- 
manifolds, by mapping the j ' = -w representation state 
onto a j = y — k state. Thus, for k = 1, in order 
to reduce j by one, one needs to initiate the system in 
the 1 000 . . .0) state (j = ^ manifold) with one injected 
bosonic excitation. After t = — Z5 7 =, the state tran- 

sitions to | j = , m = — + l), and by changing the 
relative phase between the states in the superposition 
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Total number of excitations 

FIG. 2: (Color online). Rabi splittings as a function of the 
total number of excitations n present in the TCM, for a 
nonet (j = 4) from an 8-qubit register in units of the cou- 
pling strength g (dots). We attempt at fitting these trends 
by a (8 — 2k) \fn — 2k, for k = 0, 1, 2, 3 (blue, red, green, and 
brown respectively). We see that the convergence of the eigen- 
values to the fitting trend only works for very high values of 
n > 30 and that the convergence is slowest for the highest 
eigenvalue of Hj£±. Inset zooms in on the few excitations 
regime. 

by consecutive 7V th roots of unity, a l Nl one forms the 
| j = y — l,m= — y + l) state. The application of J_ 
on that state gives 1/y/N Yli ot z N \000 . . . 0) = 0, mean- 
ing that this is the ground state of a new pseudospin 
ladder. Transformations for other values of k are more 
elaborate, with the phase shifts always carried out us- 
ing (N — k + l) th roots of unity. The details for up to 6 
qubits are presented in the supplementary material. 

Effects of decoherence. For the state-of-the-art realiza- 
tions, the resonator has coherence times exceeding these 
of the qubit by over an order of magnitude [16] . Therefore 
we explain here what happens to a three-qubit register 
if a qubit spontaneously decays. The system is initiated 
in the 1 000) <£> |1) state and evolves within the symmetric 
quartet manifold, performing Rabi oscillations with the 
frequency gy/3 between the lowest two levels of the state 
ladder. After a quarter of a period it evolves into the 
| j = §,m = -§} ® |0) = ^ (J100) + |010) + |001)) ® |0) 
state. At this point the first qubit qubits decays pbefore = 
\j = §,m = = |,m = -i| ->• pdecay resulting in 

formation of a mixed state 

Pdecay = \ | V»> (l/> | ® |0> <0| + ^ |000) (000| ® 1 0) <0| , 

IV) = |o> ® (|oi> + |io» 1 



1\ 




1 1 


to! 


2' 2/j 


2' 2 



Thus the system splits into two manifolds, one giving rise 
to the cavity line splitting (|§,— ^)|0)), and the other 
one (||,-|>|0) or (| §, + 1 §, |0)/>/2) which 
does not. Probing the transmittance of the resonator at 



its bare frequency collapses the state to either the former 
(no transmittance) with the probability of | or the latter 
(otherwise) set of states (|). In the latter case, we have 
two types of states: a dead state |000), with the proba- 
bility of Pdcad = | and a dormant one (the ground states 
of the doublets, ^dormant = f )■ In contrast to the dead 
state, the dormant one can be woken up. To do that one 
needs to perform the conditional phase shift operation 
shifting the states {ji = |,j 2 = f } -> {h = = §}• 
Upon a subsequent measurement of the cavity transmit- 
tance (at the bare frequency), negative (positive) result 
collapses the state to the j = | (j = ^) subspace. If the 
result is positive, one needs to perform one more phase 
shift operation to recover the initial state. As a result 
of this scheme one could undo the decay with the proba- 
bility of |, with no need for any direct measurements of 
the individual qubits. This process can be easily gener- 
alised to a greater size qubit register with \ j, m = —j + 1) 
states, with the recovery probability as well as the num- 
ber of potentially needed phase shift increasing with the 
increasing number of qubits. 

For initial states |000) ® \n > 1), upon a single qubit 
decay, dead and dormant states will no longer be present, 
since the ground states of the doublet can interact with 
the resonator mode, and the newly formed state with 
7i—l total excitations can form a separate quartet man- 
ifold of states. This new set of states has different Rabi 
splittings, relative to a doublet or the quartet with n exci- 
tations, allowing this more complex form of -y^n-like non- 
linearity to be used to detect any decoherence-resulting 
changes to the system as a whole. 

Strong coupling. The Jaynes-Cummings type coupling 
is only a valid description for coupling strengths of the 
order of g/£l, g/u: ~ 0.1. Greater couplings require the 
use of the counter-rotating terms in the full Rabi model 
description [17, 18]. The transition between the N spins 
state |000...0) (X) |l) boson and the state with a single exci- 
tation in the spin ensemble leads to a ^/N enhancement 
of the spin-boson coupling strength [19]. For a spin-j 
particle the single boson transition amplitude enhance- 
ment Vj, m depends on the location in the ladder, which 
ranges from y/2] — %/]V, for the states |±j) and |±j =p 1), 
to k | for \m = 0) and |m = ±1) for even, or |m = |) 
and m = — |) for odd N. Therefore, using state-of-the- 
art superconducting qubits with g/D, rj 10~ 2 [16] one 
can form a spin-10 particle which still operates at weak 
coupling [20, 21]. 

In conclusion, we have shown that a homogeneous 
Tavis-Cummings model can be simplified to produce 
higher order pseudospin interacting with a single bosonic 
mode. Such interaction, for N > 2 spins gives rise to a 
beating- like behavior. Moreover, we demonstrated that 
using a single fixed TC set-up one can emulate an en- 
tire range of pscudospins by means of conditional phase 
change transformations. We gave examples of how using 



5 




FIG. 3: (Color online). Population of the |000) state as a function of time, given that the state is initiated in the |000,n = 4} 
state. We see that the apparent beating does not have a fixed frequency. An attempt to match the beating envelope of this 
function with a cosine of a single frequency cannot be achieved for an indefinite time interval. 



such system, individual qubit decays can be detected and 
fixed. 
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Supplementary Material 
I. Details on the orthogonal block diagonalising states creation. 



A Young tableaux is a diagram composed of a collection of left-flushed boxes containing numbers which are strongly 
increasing down each column and across each row, and every consecutive row having a weakly decreasing number of 
cells (this is commonly referred to as a standard Young tableaux). By construction every such Young tableaux corre- 
sponds to an irreducible representation of a permutation group, such that this particular representation is symmetric 
under interchange of the elements in the row and anti-symmetric under the interchange of entries between columns 
(both conditions need to be satisfied simultaneously). Typically, the Young tableaux come very useful in this analysis 
and reduction of the problem into further sub-problems by means of block-diagonalisation of a Hamiltonian. 

Also here, for a homogeneous Tavis-Cummings model, with a 2 N x 2 N dimensional Hamiltonian, the very convenient 
multi-qubit states are those given by the irreducible representations of the permutation group. All tableaux that will 
be used can have at most two rows because that is the greatest number of ways in which we can anti-symmetrise the 
states. We will see that the states corresponding to different shapes of the Young tableaux will always be orthogonal, 
however states with the same number of qubit excitation corresponding to the same diagram shape, but different 
Young tableaux, will no longer be orthogonal. That is why later we will show how a deviation from the Young 
tableaux formalism further simplifies the problem. 

For 2 qubits, there are only two possible Young tableaux, one that corresponds to a triplet state, and another that 
denotes a dark singlet state. 



lfiT2l,m = 1) 



-1) 



111) 
|00) 



\m , m = 0) 
|,m = 



i 

V2 

1 
V2 



(101) 
(101) 



|io» 
|io» 



where we see that the states that the new set of states forms an orthogonal set of states. 

For three qubits there are three possible standard Young tableaux, one perfectly symmetric state | |i|2|3| ), and two 



states where the third 



w 



or a second 



qubit state is antisymmetric , while it is symmetric in the first and 



second or the first and third respectively. 



Il 


2 


3| 


, m = 


11 


2 


3| 


, m = 



= 1111) 



1_[2J3J , m 



12 3 , 7TL 



1) 



^(|011) + |101) 
^ (|001) + |010) 
1 000) 



|110)) 
|100)) 



I 


3,„ = i\ 


= 71 ( 2 I 110 >- 


|101) - 


|011)) 


~[2~ 


,m = -|\ 


= 75 ( 2 I° 01 > - 


|010) - 


|ioo)) 


i 




= 7a ( 2 I 101 >- 


|110) - 


ion)) 




,m = 


= 75 ( 2 I° 10 > - 


|001) - 


|ioo)) 



It is easy to see that any state from the [ | 1 1 2 1 3 1 ) (quadruplet j = |) manifold is orthogonal to any state in 
for x 7^ y € {2,3}, however we see that the ( |i|^J,mi ||_I|3J m2 \ = —L§ mim2 ^ indicating that the states are 



sr. 



not orthogonal. This means that the states 
alternative orthogonal set of states 



3 = 2' m 
1 

j = \,m 



I [2j 

^j^J^ and [|p^ would be coupled, and this is why we propose using 



J 





1 






1 






i 




i 

= V3 



-j= (|011) + e 2 "/ 3 |101) + e- 2 ^/ 3 |110)) , 
l y= (|100) + e 2 "/ 3 |010) + e- 2 ™/ 3 |001)) , 
(|011) + e" 2 " /3 |101) + e 2m / 3 |110)) , 
l -j= (|100) + e- 2 " /3 |010) + e 2m/3 |001)) , 



where the additional (third) label is used to distinguish the different doublet states. 
Lastly we just present the states for four qubits 

| | 1 1 2 1 3 1 4 | ,m = 2) = |1111) 



| |1|2|3|4| ,m = 1) 

| B3H3 ,m = 0) 
|1|2|3|4| , m = —1) 

frmrm , m = -2) 

, m = 



i 



12 3 



f 



i 



nm 



m = — 1 
m = 1 
m = 



I = -1 

m = 1 
m = 
i = -1 
m = 
m = 



- (|1110> + |1101) + |1011) + |0111)) 

i (|0011) + |0101) + |0110) + 1 1001) + |1010) + |1100)) 
V6 

- (|0001) + |0010) + |0100) + 1 1000)) 
|0000) 

1 



2-s/3 



(|0111) + 1 1011) + 1 1101) - 3|1110)) 



—= (|0011) + |0101) - |0110) + |1001) - |1010) - |1100>) 
V6 



1 



2^3 
1 

2^3 



(3|0001) - |0010) - |0100) - 1 1000)) 
(|0111) + 1 1011) - 3|1101) + |1110>) 



-= (|0011) - |0101) + |0110) - |1001) + |1010) - |1100>) 
V6 



1 



2V3 
1 

2V3 



(10001) - 3|0010) + 10100) + 1 1000)) 
(|0111) - 3|1011) + 1 1101) + |1110>) 



—= (-|0011) + 10101) + |0110) - |1001) - |1010) + |1100)) 
V6 



1 



2^/1 
1 

2^3 
1 

2y/3 



(|0001) + |0010) - 3|0100) + 1 1000)) 

(2|0011) + 2(1100) - |0101) - |1001) - |1010) - |0110» 

(2|1010) + 2|0101) - |0011) - |0110) - 1 1001) - |noo)) 



and show the states that corresponding to diagrams j^J b l c l and -^--^ can be orthogonaliscd to the triplets and the 
singlets 



T 1 (4) ,m = l 
T[ i \m = 
lf\m = -l 



T 2 (4) ,m = 0) = 



T 2 (4) ,m = -1 



si 



(4) 



= -(|0111) - i| 1011) - |1101) + i|1110» 



i/4 (i|0011) - |0110) + |1001) - i|1100» 



= -(— i[0001) + |0010) + i|0100) - |1000» 



= -(10111) -11011) + 11101) -11110)) 



V2 



(|0101) - |1010)) 



= 2d 0001 ) - I 0010 ) + 1 0100 ) - I 1000 )) 



= -(|0111) +i|1011) - |1101) - i|1110» 



T 3 (4) ,m = 0) = 



•• - /4 (-|0011) + i|0110) - i|1001) + |1100>) 



-3-rri/ 



(«|0001) + |0010) - i|0100) - |1000)) 
(10011) - |0110) - |1001) + |1100)) 
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1 2 3 4 5 6 7 
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0102050 14 
1 2 5 14 42 
1 3 9 28 
1 4 14 48 
1 5 20 
1 6 27 
10 7 
1 8 
1 



TABLE II: The number of distinct multiplets depending on the number of a single cavity coupled qubits. 



(4)\ _ 



2V3 



(|0011) - 2|0101) + |0110) + |1001) - 2|1010) + |noo)) 



but one could also choose the singlets to be 



S[ 4) ^ = -j= ( 

S< 4) \ = -= (|0011) +e^|0101) + e"^|0110) +e 



^ 1 1010) + 1 1100)) 
1 1001) +e^|1010) + |1100) 



We show how to construct these states in general (more than 4 qubits) in the next section. 



Switching between pseudospins 



We now show that one can define a unitary operation that allows one to switch between different types of pseudospin 
subsystems in the N qubits-resonator system. There is a general transformation which maps the i = y representation 
state onto a j = y — 1 state. Initiating the system in the 1 000 ... 0) state, the state will belong to the J = y multiplet 
(it will be the \ j, m = —j) state). Provided one excitation is present in the resonator, the state 

transfer to |j = |,m = — y + l) = l/y/N^ |000 ... 1), where the sum is over all of the N C\ ■ 
now define a transformation 



|j = f,m = -f>will 
= TV permutations. We 



T 



N 



N 



— ,m = — — 
2 ' 2 



|000 . . . 1) = 



N N 
J 2 2 



where apj is the N root of unity. An application of J_ on that state gives 1/ 



N 



000 .. . 0) 



the sum of all apf is zero. This shows that the newly obtained state is a ground state of another, j 



multiplet. In fact there will be N — 1 such transformations Tk '■ \j 

y — l,m=y — l). giving rise to N — 1 set of ground states of j 



f ,m 



= 0, since 
= £ - 1 

2 



j = y — l,m=y — J.K giving rise to iV — 1 set ot ground states ot j = y — 1. Following another injection of 
excitation one can (for sufficiently large N) again perform a transition to a state |j = y — 1, m = — y + 2) and 
attempt a similar transformation. These transformations however become slightly more complicated for N > 5. 
Alternatively, one can just perform a two-fold excitation on the i = y multiplet and perform a conditional phase 
gate sequence. For every number of qubits N this can be calculated separately, by imposing a transformation Ansatz 



-f + 1) -> l/VF£a1v|000...1) 



s 



N 



N n 

-y +2 



N(N - 1) 



^c,|00. 



11) 



such that |cj| =1 and J-\ip') = 0. One then obtains an underconstrained system of linear equations with 
5JV (N — 1) unknowns, where one of the solutions is a set of (N — 1) roots of unity. 

This process can be generalised to transformations mapping j = y) — > |j ; = y — fc), choosing a transformed state 
Ansatz of a superposition of all possible 0s and Is permutation states with equal amplitudes but different phases. 
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FIG. 4: A quantum circuit diagram, showing a manipulation of (a) the first excited quartet state to the ground state of the 
doublet state (b) the first excited quintet state to the ground state of the triplet state (c) the first excited triplet state to the 
ground state of the singlet state. Here all the gates are the phase change on 1 1} gates with a phase indicated inside the gate-box. 



Requiring that the transformed state forms the ground state of a new ladder of states, i.e. J~\j = y — £;) = 0, one 
obtains a set of N Cj~-i equations with N Ck unknowns, leaving a freedom of N Ck-\ — N Ck parameters. We can 
always set on of these parameters to 1, since the state can have an arbitrary global phase, and then we will find that 
among the multitude of possible solutions there is going to be a subset where coefficients Cj are the (N — k + l) th 
roots of unity. We present the detailed transformation and the obtained states of a six qubit register in Figure 5. 




FIG. 5: A quantum circuit diagram, showing a manipulation of the first second and the third excited septet state to the ground 
state of the quintet, triplet states or the singlet state. Here the angles 9 inside the gate boxes denote the phase shift operation. 
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II. Details of the iV-qubit propagator calculations 

Here we present two methods of calculating the j-particle propagator. We will first show the method involving the 
more purists' approach, to later present a short-cut. 

In the resonant case we can split the propagator into two classes of matrices with an even (odd) parity {Vq } ({ V\ }) , 
where a, b are just set indices. These obey Vq Vq = Vq, V" V\ = Vq, and V C L V\ = . It is very easy to find the form 
of these matrices; the even (odd) parity ones have (non) zero entries in the cells where the column and the row digit 
add up to an odd integer, and vice-versa. Put differently, the even (odd) parity matrices are the even (odd) powers 
of Hj . In what follows we will show how to solve an arbitrary qubit number model and then explicitly solve a 
three-qubit dynamical model. 

In order to find Uj (t) = exp {igtH'^ c /tij = ^ (^/ C ) we nave to calculate an arbitrary power of matrix 

/ ' wrncn i s a nontrivial task due to commutation relations of & and . This process can however be made simple 
when using the aforementioned parity decomposition. Using the simplest nontrivial even parity matrix we can then 
find a recursive relation between the even or odd matrices. Defining M = Hj for shortness, we find that 

2j+l 
P. 9=0 

We now form an Ansatz for M 2k and M 2k+1 , k € Z. In order to avoid the problems with the ordering ambiguity 
we see that the (p th , q th ) entry in the propagator must be a process involving a \p — q\ photon transfer (absorption if 
q > p and emission otherwise), moreover this amplitude depends on the number of excitations available. Therefore 
we can make the following Ansatz 

2j+l 2j 

p,q—l r—1 
2j+l 2] 

p,q—l r—0 

where /a™ (n) is an amplitude, with the labels a,p,m denoting the number of photons transfered, a function label 
(as there are many different functions with the same number of photons transferred), and the power of the matrix 
that these functions appear in. 

With this Ansatz we can now use that M 2 M 2k = M 2k M 2 = M 2 ( fc+1 ) to find the relationship between these 
functions. In this way we find that upon equating the first column of the left and the right hand side of this relation 
we get a recursive relation for the even photon transfer amplitudes in the first column of the propagator 

j2(fc+l) = J^j2k 

ftn_( Am) Am) Am) Am) \ T 

J ~ \ JOS J2,l J4,l ' ' ' J2j,l J 

3 

M = zJ S P'i ((" + ^tftp-j-i + + l + 2p)v 2 hp _ j ) + 6 p -i iq Dj i2 q-jVj,2q-j+l + 6 p , g -i (h + 2q) (n + 2q- 1) v jt 

p,q=0 

The equation above shows simply a set of linear combination of geometric series, which can be unwound by diag- 
onalising M. and determining it's eigenvalues and using the initial condition / = Yli—i Si i. These eigenvalues arc 
then the common ratios and when we sum them up the square-roots of these eigenvalues become the photon-number 
dependent multi-qubit state Rabi frequencies, however due to a number of them we see that these frequencies form 

/- \ b ' J+1 

a mutually modulated state transition pattern e.g. I Uj (t) ) = ^2 Ak cos u>kt. So the diagonalisation of this 

matrix here is what stands in a way of calculating the eigenvalues and the coefficients in front. This can be done 
analytically up to four (j = 2) identical resonantly coupled qubits where there are three frequencies, but one of them 
is zero (as it is the case every time when an even number of qubits is present). The functions present in the odd 
parity matrix can be found by identically starting with M 2 M 2k+1 = M 2k+1 M 2 = A/ 2 ( fe+1 )+ 1 , where the same set of 
frequencies necessarily needs to be found since otherwise if these non-commensurate frequencies were different in the 
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odd and the even parity sectors, then the total probability would not be conserved which is a necessary condition by 
construction of the propagator. 

In these basis the resonantly coupled interaction picture Hamiltonian reads 



#j=3/2 



/o 

V 





and as a consequence 



H 



j=3/2 



= M = 



\/3a 





\ 





2 a 





2a) 





•\/3a 







o / 







2\/iaa 


7n 


+ 4 






/3n + 3 

7h + 4 2y/3aa 

2y/2,a)a) 7n + 3 

\0 2 % /3a t at 3n 

Let's define an Ansatz for the even parity matrix element 



M 



2 A- 



/ fH 





f (m) - ~ 
J2,l aa 








f (m) 
^0, 2 





j (m) « « 
J2,2 a « 







f (m) 
JO, 3 





\o 







f (m) 
JO, 4 



Then using that M 2k M 2 = M 2 ( fe+1 ) we get that 

f$ +1) = S(n+l)f^+2V3(n + l)(n + 2)f^ 
/ 2 ( X +1) = 2>/3/#> + (17 + 7n)/£> 
/ should explain how to get this result 

(m) _ A!"(A_ - A + - 4 (n + 2) - 6) + A"(A_ - A + + 4 (n + 2) + 6) 



Jo,i 

(rn) _ 2V3(A m -A 

J2,l — 



2(A_ - A+^ 



A_ - A + 

A± = 5 (n + 2) ± \J 16 (n + 2) 2 + 9 
There is also an equation which feeds on solutions for f^™ from the above set of equations 

/oS +1) (n) - (7n + 3) /£> (») + 2^3^/^ (n) aa 



(2) 



(3) 



(4) 



(5) 



where we reinstall the explicit photon number dependence, because now using g (n) a = ag(n— 1), we can rewrite 
the equation above 



with a solution 



/o5 +1) W = (7n + 3) /ff (n) + 2^ (n - 1) n/#> („ - 2) 
(m) _ <p m (4n + y_ - <p + + 6) - yy(4n - y_ + y+ + 6) 

Jo, 3 W — 



2(</>- - </>+) 
(,5-1- = 5n ± \/ 16n 2 + 9 

Similarly we have another set of equations 

ft2 +1) = (7n + 4)/^+2V3(n + l)(n + 2)/iJ 
/ 2 ( ™ +1) = 2>/3/<?>+3(n + 2)/<J 
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with solutions 



Am) 

JO, 2 



Am) 
J2,2 



V m (n- - fx + + 4 (n + 1) - 6) + m+ (/J- - M+ - 4 (n + 1) + 6) 



2(m- - 



2V3( 



and then the finally 



with a solution 



: 5(n+l)±yi6(n+l) z + 9 

/<^ +1) in) = W$ (n) + 2^^/^ (ft) aa 

= 3nf$(n)+2V3(n- 1) »/£> (n - 2) 

«;"(«;_ - k + - 4 (n - 1) + 6) + - k+ + 4 (n - 1) - 6) 

2(k_ — K + ) 



/ ( ? W 

k± = 5(n- 1)± ^16 (n- l) z + 9 
Similarly an the odd subspace is obtained, by first starting with the Ansatz 



M 2k+1 



( 




































o / 



with solutions 



f (m) 



Am) 
J3,l 



Am) 
Jl,2 



Am) 
J1.3 



V3 (A'"(A_ - A + + Ap - 6) + Ay (A_ - A+ - 4p + 6)) 
2(7p-3)(A_-A+) 
6 (A™ - Ay) 6 (A™ - Ay) 
(7p-3)(A_-A+)(7p-3)(A_-A + ) 

2 (m™ - /uy) 

M- - M+ 

y^(4g - y_ + ip+ - 6) + y>y (-4g - g_ + y>+ + 6) 
2V3(g + l)(^ + -^_) 



(6) 



In all of the expressions above, upon summing, every m th power of A± , ip± , fi± or re± gets replaced by a sine or a 
cosine of a square-root of this quantity (with some minor changes with the two and three photon entries, since the 
lowest entries there are missing when summing up). 



Fo,i 
Fo,2 

Fo,3 
Fq,4 
F2,l 



(A_ 2 + 4 (h + 2) - A + ,2 + 6) cos (u+, 2 t) - (-A_ 2 + 4 (n + 2) + A +;2 + 6) cos (u-, 2 t) 

2(A_, 2 - A +>2 ) 

(A_,i - A+,i + 4 (n + 1) - 6) cos (u)-,it) + (A_,i - A+,i - 4 (n + 1) + 6) cos (u)+,it) 

2(A_,i-A+,i) 

(A_ j0 + 4n - A + , + 6) cos (w_,o*) - (— A_,o + 4n + A +;0 + 6) cos (w +i o*) 

2(A_,o — A+,o) 

(-4 (n - 1) + A_ _i - A+ _i + 6) cos (w_ _it) + (4 (ra - 1) + A_ _i - A+ _ x - 6) cos (w+ -it) 

2(A_,_i - A+,_i) 

» (\/A-,2(A-, 2 - 4 (n + 2) - A + , 2 - 6) sin (u-, 2 t) + y/A7^(A_, 2 + 4 (n + 2) - A + , 2 + 6) siu ju+tf)) 

2V3(p-l)(A_, 2 -A +i2 ) 
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Fl,2 
-Pi, 3 
-^2,1 
^2,2 
■^3,1 



2i (*/A_ ( i sin (ui-jt) — V^+,i sin (w+^f)) 
A-,i — A + x 

z ( v /A_,o(A_,q - An - A +;0 + 6) sin (u;-,of) + t/A+^A-.q + An - \ +:0 - 6) sin (u + , t)) 

2V3(n+l)(A_, - A+, ) 

2y3 (cos (w-^i) — cos (w+^t)) 

A-,2 — A+,2 
2\/3 (cos — cos (w +! it)) 

A-,i — A+_i 

j ( yXr^(-A_, 2 + 10 (n + 2) + A + , 2 ) sin Qj_ 2 *) + V /A~^(-A_^ - 10 (ft + 2) + A + , 2 ) sin (u + , 2 t)) 

3(p 2 -l)(A_, 2 -A + , 2 ) 



9y/^±,x ^±,x — 5 (n + sc) ± \Jl6 (n + x) 2 + 9 



(7) 



A short-cut method uses the fact that the system conserves the total number of excitations. Assuming, for now, 
that we start in a state \j, m = j, n) we see that there is only a finite number of states that will be populated, and 
every consecutive state down the j— ladder will have one more photon. We can then write (using the three qubit 
example) the Hamiltonian as 



H 



j=3/2 



\ 



( Via 

Via) 2a 

2a) Via 

\0 Via) J 



/o 

VlVn~- 


Vo 



VIVr 



1 



2Vn 




+ 1 

2Vn~ 
~2 



VlVn + 2 
Vlyfn + Z j 



(8) 



and then the propagator simply takes the form exp igtH'j_ 3 ; 2 > which can be readily calculated using available com- 
puter algebra software. Upon obtaining the result one has to match the expressions found from matrix exponentiation, 
to the structure of the propagator Ansatz, i.e. identify the entries in the matrix in the pseudospin space where the 
creation or an annihilation operator should be placed. Afterwards the propagator needs to be subdivided into regions 
shown below 



/ a a a a\ 
a b b b 
abac 

\ a b c d ) 



(9) 



such that the pattern continues for larger matrices introducing regions e, /, etc. Subsequently, in the first row and 
first column (denoted with a), every instance of Vn + k (k being an integer) in a numerator should be treated as a 
creation or annihilation operator, which one needs to place in a normal ordering convention i.e. write the annihilation 
(creation) operators to the right (left) of the number operators. Similarly, in the second region (denoted with b), one 
performs the same process, however after placing the creation and annihilation operators, every remaining instance 
of variable n needs to be replaced with n — 1. Then for region c the replacement is n — > n — 2, for region d it is 
n — > n — 3 etc. 



III. Solutions to the system with dephasing 

The same method as outlined above can be used when studying the solutions to the master equation of a resonator 
coupled to a qubit which undergoes dephasing. The equation of such a system in the interaction picture reads 

p = -i [Sa z + g (o + a + u_a t ) ,p] + - [a z pa z - p) . (10) 

The perceived difficulty of this problem is two-fold. For one, the system is highly nonlinear, where the number of 
bosonic operators in a nonlinear fashion affects the transition amplitude of the qubit state (also before we saw h 
present in a square-root as an argument of a sine or a cosine). Secondly, the evolution is no longer unitary and it 
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FIG. 6: (Color online). Population of the 1 000} state as a function of time, given that the state is initiated in the |000, n = k) 
state, where k — to k = 8 increasing left to right and down the column. 



requires density operator formalism treatment and a use of superoperators acting on the density matrix from both 
sides simultaneously. 

We can address the first problem by taking the density matrix in a general form p = ^2\i, n)(j, m\, where i (n) and 
j (m) are the qubit (resonator) states, and rewriting it in terms of a product state p = Y1 K n )|j\ TO )- As a result of 
that for any A\pB 2 — >• A\ <3> B 2 p, where we have adapted a notation where the subscripts denote the space that the 
operators arc acting on 

ai\i,n)\j,m) = y/ri\i, n - l)\j, m) 
a 2 \i,n)\j,m) = \Jm + n)\j 7 m + 1) 
a\\i, n)\j, m) = \fn +T|i, n+ l)\j, m) 
a\\i,n)\j,m) = y/m\i, n) \j, m - 1) 

With this formalism in mind the equation of motion above takes the form 

iga 2 



iga\ 



Hp 







-iga\ 



—igai 
iS 

—(p — iS 
iga 2 



-igai 
iga\ 





(Pn\ 




Pw 




Poi 


) 


K Poo J 



(11) 



with pij denoting the qubit space indices of the density matrix, leaving the photonic ones implicit. The general form 
solution to the equation above reads p (t) = exp (jltj p (0), which again can be found in a closed form provided that 

we will be able to find H k for any k. For the time being we will analyse the case of 5 = 0. 
By analysing the first few powers of H and using the commutation relations 



a] , a. 



M-ir 



(12) 
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(13) 



we see that the matrices H k will always have a particular structure with allowing us to form an Ansatz 

f fi (&) fs (k) 4 fs ( k ) "i ho ( k ) «i4 \ 

jjk _ a 2 f 5 (k) h a 2 f 6 (k) oi f 9 (k) oi 

4/s (fc) 4/ 6 (fc) 4 / 3 (fc) / 7 (fc)4 

\ O1S2/10 (fc) 4/9 (k) a 2 f 7 (k) f 4 (k) J 
where functions fj are functions of the number of bosonic excitations n, m in either of the spaces. This allows us 
to form recursive equations based on H k+1 = HH k , which are linearly coupled geometric progressions which can be 
uncoupled by a simple diagonalisation process. 

°° (if) k 

The solutions fj to these recursive equations, when summed up to obtain Fj = W form a product space 

k=0 

propagator with the same form as the Ansatz above, but with fj — > Fj, where Fj read 



Fi = 

F 2 = 

F 3 = 

F A = 

F 5 = 

F e = 

F 7 = 

Fs = 

F 9 = 

F10 



X 



/± ,x,y 



-e-T(c_ 1,1 +C+.i,i + 



1 _t4> 

r~ 2 
1 

r~ 2 
1 t<t> , 

-e 2 (C_,o,o + C+,0,0 + ^-,0,0 + £+,0,0) 



»g _M / 1,2 



(C-,0,1 + C+,0,1 ^ ^-,0,1 — S+,0,1) 

(C-,1,0 + C+,1,0 — £-,1,0 — S+,l,o) 



e"* (C-,1,1 - C+,1,1 - + 

2^/(^1 + l)(n 2 + l) 
*5 -*£ ( \,1 Q . 1,2 c 



1,0 



«3 _M f 2,1 
-— e 



r?-,!,!^-,!,! + 7 7+,l,l^+-l,l 
«ff -if. f 2,1 Q . 2.1 c \ 

~~0 e [ t I-,o,i s -,oa + V+,o,iS+,o,iJ 

e~~%~ (C-,1,1 - C+,1,1 + ff-,1,1 - 5+, 1,1) 

2^/(^1 + l)(n 2 + l) 
0sin (lfftA+,3^) 



C±,x,i/ = cos ( -gt\ ±tX , y 



1± 



"9 + 2/ 

In the < — > 00 limit, the propagator above takes a form 

1 



}j^ exp \ Ht i = o (IphXphI + IpooXpooI) 



4«2 



2 V y/fc + 1) (n 2 + 1) 



110 



Assuming that we initialise the cavity-qubit system in a pure state pi n i t = \ip)(ip\, where = cos 011, ft) + 
e la sin#|+, m), then under the evolution of this equation of motion the state will become 

Pfinai = -diag{c|n)(n| + s|m + l)(m + 1|, 
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c\n + l)(n + 1| + s|m)(m|} 

where c = cos 2 9 and s = sin 2 8. This means that regardless of what the initial state of the qubit-cavity system is, 
the qubit dephasing will lead to loss of any coherence or entanglement within and across the subsystems. Also when 
we trace out the cavity the qubit is left in the state pa na i = \ (|t)(tl + I4)(-J-I)i w hich is different if the coupling to the 
resonator was absent cos 2 0|t)(tl + sm2 (4-1- 

If we treat the detuned case, then based on the same Ansatz and diagonalising the recursive relation, we need to 
find the eigenvalues (which are the common ratios of the geometric progression) of the coupled system of functions 
from the first column of the Ansatz. These eigenvalues later act as the modified Rabi frequencies, and these can be 
found by means of finding the roots of the characteristic polynomial 

P (x) = x 4 + 2cj)x z + (2n ig 2 + 2n 2 g 2 + 4 5 2 + S 2 + 4> 2 ) x 2 

+ (40 5 2 + 24>n ig 2 + 2§n 2 g 2 ) x(j) (14) 

which in the limit of 5 — > has the roots —cf> + i\±,i t i and its complex conjugate. For non-zero d, to the best of the 
authors knowledge, there arc no analytical solutions. 



